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Abstract. 

In this paper we consider from the point of view of differential geometry and of the 
theory of integrable systems the so-caUed WDVV equations as defining relations of 2- 
dimensional topological field theory. A complete classification of massive topological con- 
formal field theories (TCFT) is obtained in terms of monodromy data of an auxiliary 
linear operator with rational coefficients. Procedure of coupling of a TCFT to topologi- 
cal gravity is described (at tree level) via certain integrable bihamiltonian hierarchies of 
hydrodynamic type and their r-functions. A possible role of bihamiltonian formalism in 
calculation of high genus corrections is discussed. As a biproduct of this discussion new 
examples of infinite dimensional Virasoro-type Lie algebras and their nonlinear analogues 
are constructed. As an algebro-geometrical applications it is shown that WDVV is just the 
universal system of integrable differential equations (high order analogue of the Painleve- 
VI) specifying periods of Abelian differentials on Riemann surfaces as functions on moduli 
of these surfaces. 
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This paper is an extended version of the talk given at the Workshop "Integrable 
Systems" (Luminy, July 1991). Also more recent results [39] are included. I dedicate it to 
the memory of J.-L.Verdier. 

Introduction. 

A quantum field theory (QFT) on a D-dimensional manifold M consists of: 

1) . a family of local fields ^^(x), x E M (functions or sections of a fiber bundle over 
M). A metric gij{x) on M usualy is one of the fields (the gravity field). 

2) . A Lagrangian L = L{(p, (p^, ...). Classical field theory is determined by the Euler 
- Lagrange equations 

3) . Procedure of quantization usualy is based on construction of an appropriate path 
integration measure [d(/)]. The partition function is a result of the path integration over 
the space of all fields (f){x) 




Correlation functions (non normalized) are defined by a similar path integral 

< Mx)My) ■ ■ ■ >M= J [d4>]Mx)My) ■ ■■e'^^'^l 

Since the path integration measure is almost never well-defined (and also taking in account 
that different Lagrangians could give equivalent QFT) an old idea of QFT is to construct 
a self-consistent QFT by solving a system of differential equations for correlation func- 
tions. These equations were scrutinized in 2D conformal field theories where D=2 and 
Lagrangians are invariant with respect to conformal transformations 

6gij{x) = egij{x), 6S = 0. 

This theory is still far from being completed since complexity (and, probably, nonintegra- 
bility) of the differential equations determining correlators. 

Here I will consider another class of solvable QFT: topological field theories. These 
theories admit topological invariance: they are invariant with respect to arbitrary change 
of the metric gij{x) on M 

dgij{x) = arbitrary, dS = 0. 

On the quantum level that means that the partition function Zm depends only on topology 
of M. All the correlation functions also are topological creatures: they depend only on the 
labels of operators and on topology of M but not on positions of operators 

< 4'a{x)(f)f3{y) . . . >M = < (t>a4>f3 ' ' ' >M ■ 

The simplest example is 2D gravity with the Hilbert - Einstein action 




= Euler characteristic of M. 
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There are two ways of quantization of this functional. The first one is based on an ap- 
propriate discrete version of the model (M — > polihedron). This way leads to considering 
matrix integrals of the form [55] 



(t) = / exp{-tr(X2 



Jx*=x 



+ tiX^ + t2^^ + . . .)dX 



where the integral should be taken over the space of all N x N Hermitian matrices X. 
Here ti, t2 ■■■ are called coupling constants. A solution of 2D gravity [1] is based on the 
observation that after an appropriate limiting procedure ^ oo (and a renormalization 
of t) the limiting partition function coincides with r-function of the KdV-hierarchy. 

Another approach to 2D gravity is based on an appropriate supersymmetric extension 
of the Hilbert - Einstein Lagrangian [2]. This reduces the path integral over the space 
of all metrics gij{x) on a surface M of the given genus g to an integral over the finite- 
dimensional space of conformal classes of these metrics, i.e. over the moduli space Mg 
of Riemann surfaces of genus g. Correlation functions of the model are expressed via 
intersection numbers of some cycles on the moduli space [2-4, 46, 49] 



(here the supscript g means correlators on a surface of genus ^r). This approach is often 
called cohomological field theory. 

It was conjectured by Witten that the both approaches to 2D quantum gravity should 
give the same results. This conjecture was proved by Kontsevich [42-43] . He showed that 
the generating function 



(the free energy of 2D gravity) is logarythm of r-function of a solution of the KdV hierarchy 
(this was the original form of the Wittcn's conjecture). The r-function is specified by the 
string equation (see eq. (3.16b) below). 

Other examples of 2D TFT constructed in [2-6, 8-9, 46-49, 57] proved that these 
could have important mathematical applications, probably being the best tool for treat- 
ing sophisticated topological objects. In these examples correlators can be expressed via 
intersection numbers on moduli spaces (or their coverings [48]) of holomorphyc maps of 
Riemann surfaces to a complex (or even almost complex) variety { topological sigma-models 
[2]) or via intersection form of a singularity in the catastrophe theory {topological Landau - 
Ginshurg models [9, 7]; see also [10]). (We do not discuss here interesting relations between 
these models.) This gives rise to the following 

Problems. What could be an intrinsic origin of integrability in 2D TFT? How one 
can classify 2D TFT? Is it possible to find an analogue of the KdV hierarchy for calculating 
the partition function of a given TFT model? 



(I)a ^ Ca e H^{Mg), a e N 



< 0a0/3 • • • >s= #(Ca n n . . .) 




3 



In this paper an approach to these problems is proposed being based on differential 

geometry and on the theory of classical integrable systems of KdV type. Main ingredient 
of my approach is Hamiltonian formalism of integrable hierarchies of KdV type (see, e.g., 
[21, 25, 29]) and, especially, Hamiltonian analysis of semi-classical limits of these systems 
[18-21]. 

Let me start with considering matter sector of a 2D topological field theory. That 
means that the set of local fields (pi{x) , . . . , (pnix) (the so-called primary fields of the 
model) does not contain the metric. (Afterwards one should integrate over the space of 
metrics. This should give rise to a procedure of coupling to topological gravity that will 
be described below.) Then the correlators of the fields 0i(a;), (j)n{x) obey very simple 
algebraic axioms (a consequence [51] of the general Atiyah's axioms [50] of a topological 
field theory). 

Let 

VaP =< <Pa4>/3 >0 

(0 means genus zero correlator), 

Ca/37 =< (f>a(l>l3(f>'r >0 

Then 

1) These tensors are symmetric and det{r]ai3) ^ 0. I will use the tensor r]ai3 and the 
inverse (r/"^) — (rya/j)""*^ for lowering and raising indices. 

2) ^ap ~ V'^^^aiie is a tensor of structure constants of a commutative associative algebra 
A with a unity. That means that for a basis ei, e„ in A the multiplication law has the 
form 

7 

(We will normalise a basis in such a way that ei = the unity of A. So c^^, = 6^.) 

3) Let H — rj^^Caep G A. Then for correlators of genus g the following formula holds 

< (j)a ■ ■ - (j)-^ >g = < e^. . .e^,H^ > . 

On this way 

Topologicaly invariant Lagrangian — > correlators of local physical fields 

we lose too much relevant information. To capture more information on a topological 
Lagrangian we will consider a topological field theory together with its deformations pre- 
serving topological invariance 

(t" are coupling constants). Here we use ideas and results of [8, 51]. In these papers it was 
proposed a general construction of a class of 2D TFT by twisting of N=2 superconformal 
field theories. So-called topological conformal field theories (TCFT) are obtained by this 
procedure. For any TCFT with n local observables (primary operators) it was constructed 
a canonical n-parameter deformation preserving topological invariance. All the correlators 
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of the primary fields (pi, 0^ in the perturbed TCFT now depend on couphng parameters 
ti, tn- This dependence is not arbitrary but obeys the following eguations: 

1) . r]ocf3 = const in t 

2) . cff,^5^ 

3) . Cap-y — qILq-ij3q-i--i for some function F{t) (primary free energy). 
Equations of associativity give a system of nonlinear PDE for Fit) 

d'Fjt) d^Fjt) ^ d^Fjt) d^Fjt) 

with the constraint 

^ ' Vap. (0.2) 



dt^dt^dtf^ 

These equations were called in [39] Witten - Dijkgraaf - E. Verlinde - H. Verlinde (WDVV) 
equations. In fact in TCFT one should assume invariance of a solution with respect to 
scaling transformations of the form 

VaP ^ c'^-^'^'SaP 
CaP^C^-^''-''-clp 

for some numbers [charges of the fields (/)q) and d [dimension of the model). 
My program now is: 

1. To classify 2D TFT as solutions of WDVV equations, and 

2. For any solution of WDVV (I recall that this describes the matter sector of a TFT 
model) to construct (i.e., to calculate partition function and correlators) a complete TFT 
model (coupling of the given matter sector to topological gravity). 

The problem 1 was investigated in [39] . A Lax pair for the WDVV equations was con- 
structed. For the so-called massive TCFT models where the Frobenius algebra for almost 
all t has no nilpotents it was shown that solutions of the WDVV equations form a "-('^^~^) _[_]__ 
dimensional family (where n is the number of primaries). It turns out that the WDVV 
equations for massive case are equivalent to equations of isomonodromy deformations of 
an ordinary linear differential operator with rational coefficients. These isomonodromy 
deformation equations coincide with the Painleve-VI equation (for n = 3) and with high 
order analogues of the Painleve-VI for n > 3. Monodromy data (i.e. Stokes matrices) of 
the operator with rational coefficients serve as parameters of massive TCFT-models. 

Concerning the second problem my conjecture is that the set of solutions of WDVV 
parametrizes a big class of hierarchies of l-|-l-integrable systems. All the wellknown hier- 
archies are in this class but they are only isolated points in it. 

The basic idea of construction of these integrable hierarchies comes from the standard 
in quantum field theory Feynmann diagram expansion machinery. In 2D QFT it becomes a 
representation of partition function and correlators as a sum of contributions from surfaces 
of different genera g [genus expansion). The idea is that the genus expansion of the 
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partition function should coincide with the small dispersion expansion (see below) of the 
T-function of some integrable hierarchy. 

A first step on this way has been done in [39] : for any solution of WDVV a hierar- 
chy of integrable Hamiltonian equations of hydrodynamic type was constructed such that 
r-function of a particular solution of this hierarchy coincides with the genus zero approx- 
imation of the correspondent TFT model coupled to gravity (see Section 3 below). From 
the point of view of WDVV equations the hierarchy determines a family of symmetries 
of the equation (0.1) (see below, Proposition 3.1). To go further one should solve a non- 
standard "inverse problem": to reconstruct integrable hierarchy from the zero-dispersion 
limit of it. Some examples of such a reconstruction are discussed in Sections 3, 4. Probably, 
bihamiltonian formalism could be useful to complete solution of this problem. 

Examples of solutions of WDVV and of corresponding integrable hierarchies are de- 
scribed in Section 4. Almost all the known examples are obtained as a result of analysis of 
the semiclassical (particularly, dispersionless) limiting procedure in integrable hierarchies 
of KdV type. More precisely, let 

dt^y" = fkiy, d^y, dly, . . .), a = 1, . . . , /, A; = 0, 1, . . . 

be a commutative hierarchy of Hamiltonian integrable systems of the KdV type. "Hierar- 
chy" means that the systems are ordered, say, by action of a recursion operator. Number of 
recursions determine a level of a system in the hierarchy. Systems of the level zero form a 
primary part of the hierarchy (these correspond to the primary operators in TFT); others 
can be obtained from the primaries by recursions. The hierarchy posesses a rich family of 
finite-dimensional invariant manifolds. Some of them can be found in a straightforward 
way; one needs to apply sophisticated algebraic geometry methods [28] to construct more 
wide class of invariant manifolds. Any of these manifolds after an extension to complex 
domain turns out to be fibered over some base M (a complex manifold of some dimension 
n) with m-dimensional tori as the fibers (common invariant tori of the hierarchy). For 
m = M is nothing but the family of common stationary points of the hierarchy. For 
m > M is a moduli space of Riemann surfaces of some genus g with certain additional 
structures: marked points, marked meromorphic function etc. These are the families of 
finite-gap ("(/-gap") solutions of the hierarchy. The main observation is that any such M 
determines a solution of WDVV equation (M is a Frobenius manifold in the terminology 
of Section 1 below, or the "small phase space" of a TFT theory in the terminology of [3-4, 
46]). For m = and the set M of stationary points of the Gelfand - Dickey hierarchy this 
essentialy follows from [8, 11]; for general case (including arbitrary genera g) a construction 
of solution of WDVV was given in [13-14] (see also recent preprint [44]). 

To give an idea how an integrable Hamiltonian hierarchy of the above form induces 
tensors cj^^, rjafs on a finite dimensional invariant manifold M I need to introduce the 
notion of semiclassical limit of a hierarchy near a family M of invariant tori (sometimes it 
is called also a dispersionless limit or Whitham averaging of the hierarchy; see details in 
[15-21]). In the simplest case of the family of stationary solutions the semiclassical limit 
is defined as follows: one should substitute in the equations of the hierarchy 

X ^ ex = X, tk ^ etk = Tk 
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and tend e to zero. For more general M (family of invariant tori) one should add averaging 
over the tori. As a result one obtains a new integrable Hamiltonian hierarchy where the 
dependent variables are coordinates , v"' on M and the independent variables are the 
slow variables X and Tq, Ti, ... . This new hierarchy always has a form of a quasilinear 
system of PDE of the first order 

dT,vP^CkPgiv)dxv'i, A; = 0,1,... 

for some matrices of coefficients Ck^{v). One can keep in mind the simplest example of 
a semiclassical limit (just the dispersionless limit) of the KdV hierarchy. Here M is the 
one-dimensional family of constant solutions of the KdV hierarchy. For example, rescaling 
the KdV one obtains 

ut = uux + e^uxxx 
(KdV with small dispersion). After e — > one obtains 



Ut = uux- 



The semiclassical limit of all the KdV hierarchy has the form 

dTkU = -TT^xu, A; = 0, 1, 

A semiclassical limit of spatialy discretized hierarchies (like Toda system) is obtained 
by a similar way. It still is a system of quasilinear PDE of the first order. 

Let us come back to determination of tensors rya/j, c^^^ on M. To introduce ?7a/3 
we need to use the Hamiltonian structure of the original hierarchy. A semiclassical limit 
(or "averaging") of this Hamiltonian structure in the sense of general construction of 
S.P.Novikov and the author induces a Hamiltonian structure of the semiclassical hierarchy: 
a Poisson bracket of the form 

{v^{X), ^^(y)}semiciassicai = g^\v{X))[5ldx5{X -Y) - Vl,{v)v'^5{X - Y)] 

where g^'^{v) are contravariant components of a metric on M and r|^(f ) are the Christoffel 
symbols of the Levi-Civita connection for g^'^{v) (the so-called Poisson brackets of hydro- 
dynamic type) . (Strictly speaking the metric and the connection are defined on a real part 
of M that parametrizes smooth solutions of the original hierarchy with some reality con- 
straints. But the formulae for the metric and the connection admit an extension onto all 
M.) From the general theory of Poisson brackets of hydrodynamic type [18-21] one con- 
cludes that the metric g^'^{v) on M should have zero curvature. So local fiat coordinates 
t^, on M exist such that the metric in this coordinates is constant 

dt" dt^ , 

d^d^,9'"'{v) = V.p = const. 
The Poisson bracket { , } semiclassical in these coordinates has the form 

{r(X),t^(y)}semiclassical = r/"^(5'(X - Y). 
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The tensor (r/a/j) = {r]°''^)~^ together with the flat coordinates is the first part of a 
structure we want to construct. (The flat coordinates t^, can be expressed via 

Casimirs of the original Poisson bracket and action variables and wave numbers along the 
invariant tori - see details in [18-21].) 

To define a tensor c^^(t) on M (or, equivalently, the "primary free energy" F{t)) we 
need to use a semiclassical limit of the r-function of the original hierarchy [11, 53-54, 61] 

logrsemiciassicai(7b,Ti, . . .) = lime"^ logT(eto, eti, . . .). 

e— >0 

Then 

F — ^semiclassical 

where Tsemiciassicai should bc cousidcrcd as a function only on n of the slow variables of 
the same level. (To satisfy the normalization (0.2) one should choose properly these n 
slow variables and normalize values of others. This specifies uniquely the semiclassical 
T-function.) The semiclassical r-function as the function of all slow variables coincides 
with the tree-level partition function of the matter sector rjafs, c^p coupled to topological 
gravity. 

Summarizing, we can say that a structure of Frobenius manifold (i.e., a solution of 
WDVV) on an invariant manifold M of an integrable Hamiltonian hierarchy is induced by 
a semiclassical limit of the Poisson bracket of the hierarchy and of the r-function of the 
hierarchy. So the above conjecture can be reformulated as follows: WDVV equations just 
specify the semiclassical limits of r- functions of Hamiltonian integrable hierarchies. 

I do not consider in this paper one more type of integrable systems being involved 
in 2D TFT: the so-called equations of topological- antitopological fusion proposed in [40]. 
These equations describe the ground state metric on a given 2D TFT model. See [45] about 
the theory of integrability of these equations. An interesting relation of these equations to 
the theory of harmonic maps also was found in [45] . 

1. Geometry of Frobenius manifolds. 

I recall that A is called a Frobenius algebra (over R or C) if it is a commutative 
associative algebra with a unity and with a nondegenerate invariant inner product 

< ab,c >=< a,bc> . (1.1) 

If e is the unity of A then the invariant inner product on A can be written in the form 

< a,b >— u!e{ab) (1.2a) 

where 

u!e{a) =< e,a> . (1.26) 
Moreover, for any linear functional u E A* the inner product 

< a,b >u!= cij{ab) (1.3) 
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is invariant. It is nondegenerate for generic u. Any invariant inner product on a finite- 
dimensional Frobenius algebra A (only finite-dimensional algebras will be considered) can 
be represented in the form (1.3). 

If Ci, i = 1, . . . , n is a basis in A then the structure of Frobenius algebra is specified 



by the coefficients rjij , c^- where 



< Ci^Cj >= rjij (1.4a) 



eiCj = c'y.ek (1.46) 



(summation over repeated indices will be assumed). The matrix rjij and the structure 
constants c^^ satisfy the following conditions: 

r]ji = r]ij, det{r]ij)^0 (1.5a) 

^ij'^ik — ^is'^jk (1.56) 

(associativity) , 

(^ijk — Vis^^jk — ^jik C^fcj (1.5c) 

(commutativity and invariance of the inner product). If e = (e*) is the unity of A then 

e^cl^ = S] (1.5d) 

(the Kronecker delta). 

1-dimensional Frobenius algebras are parametrized by 1 number (length of the unity). 
Any semisimple n-dimensional Frobenius algebra is isomorphic to the direct sum of n 
one-dimensional Frobenius algebras 

fifj = Sijfi, < fi, fj >= r}ii5ij. (1.6) 

Moreover, any Frobenius algebra without nilpotents is a semisimple one. 

Let us consider a particular class of deformations of Frobenius algebras. 

Definition 1.1. A manifold M is called quasi- Frobenius if it is equipped with three 
tensors c = (c^j(x)), ij = {riij{x)), e = {e^{x)) satisfying (1.5) for any x G M. 

In other words these three tensors provide a structure of Frobenius algebra in the 
space of smooth vector fields Vect{M) over the ring T{M) of smooth functions on M: 

[v ■ w]^{x) = c'lj{x)v\x)w^{x), (1.7a) 
<v,w> {x) = r]ij{x)v^{x)w-' (x) (1-76) 

for any v, w E Vect{M). 

Complex quasi-Frobenius manifolds are defined in a similar way but the tensors c, r], 
e should be holomorphic. They provide a structure of Frobenius algebra in the space of 
holomorphic vector fields over the ring of holomorphic functions. 

Informaly speaking, n-dimensional quasi-Frobenius manifolds are n-parameter defor- 
mations of n-dimensional Frobenius algebras. For any x E M the tangent space Ta;M is a 
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Frobenius algebra with the structure constants c^j{x), invariant inner product 'r]ij{x), and 
unity e*(a;). 

As it was explained above, in physical applications there are additional restrictions 
for quasi-Frobenius manifolds. 

Definition 1.2. A quasi-Frobenius M is called Frobenius manifold if the invariant 
metric 

ds"^ = r]ij{x)dx'^dx^ (1.8a) 
is flat, the unity vector field e is covariantly constant 

Ve = (1.86) 

(here V is the Levi-Civita connection for ds'^) and the tensor 

Vz < u ■ v,w > (1.8c) 

is symmetric in the vectors u, z. 

Localy Frobenius manifolds are in 1-1 correspondence with solutions of WDVV equa- 
tions (i.e., with 2D TFTs). Indeed, for the flat metric (1.8a) localy flat coordinates 
exist such that the metric is constant in these coordinates, ds"^ — rjotpdt^'dtl^ , rjap — const. 
The covariantly constant vector field e in the flat coordinates has constant components; 
using a linear change of the coordinates one can obtain e'* = 5f . The tensor CaP'^it) in 
these coordinates satisfles the condition 



dsCapj = (1.9a) 
This means that CQ,/3^(t) can be represented in the form 

cap^it) = d^d(3d^F{t) (1.96) 

for some function F{t) satisfying the WDVV equations. 

The flrst step in solving WDVV is to obtain a "Lax pair" for these equations. The most 
convenient way is to represent them as the compatibility conditions of an overdetermined 
linear system depending on a spectral parameter A. 

Proposition 1.1. A quasi- Frobenius manifold is Frobenius iff the unity e is covari- 
antly constant and the pencil of connections 

Vu{>)v = VuV^-Xu-v, u,veVect{M) (1.10) 

is fiat identicaly in A. 

Flatness of the pencil of connections (1.10) is equivalent to flatness of the metric rj 
and to the equation 

Vu{v ■ w) — V v{u ■ w) + U ■ WvU! — V ■ Vulf = [u, v] ■ w (l-H) 

for any three vector flelds u, v, w. Here [u, v] means the commutator of the vector flelds. 
This equation is equivalent to the symmetry of the tensor (1.8c). 
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Corollary. WDVV is an integrable system. 

Indeed, WDVV is equivalent to compatibility of the following linear system 



Va(A)C = 0, Q! = l,...,n, 



(1.12a) 



(here ^ is a covector field), or, equivalently, in the flat coordinates t 



a 



(1.126) 



Compatibility of the system (1.12) (identicaly in the spectral parameter A) together with 
the symmetry of the tensor Ca/s-y — fJaeCpj is equivalent to WDVV. 

It turns out that symmetries of Frobenius manifolds play an important role in geomet- 
rical foundation of TFT. We start with the notion of algebraic symmetry of a Frobenius 
manifold. 

Definition 1.3. A diffeomorphism / : M — > M of a Frobenius manifold is called 
algebraic symmetry if it preserves the multiplication law of vector fields: 



(here /* is the induced linear map /* : T^M — Tf(^x)M). 

Proposition 1.2. Algebraic symmetries of a Frobenius manifold form a finite-dimen- 
sional Lie group G{M). 

The generators of action of G{M) on M (i.e. the representation of the Lie algebra of 
G{M) in the Lie algebra of vector fields on M) are the vector fields w such that 



for any vector fields tt, 

Note that the group G{M) always is nontrivial: it contains the one-parameter sub- 
group of shifts along the coordinate t^ . The generator of this subgroup coincides with the 
unity vector field e. 

The group G{M) can be calculated for the important class of massive Frobenius 
manifolds. 

Definition 1.4. A Frobenius manifold is called massive if the algebra on T^M is 
semisimple for any x E M. 

In physical language massive Frobenius manifolds are coupling spaces of massive TFT 
models. 

Main lemma. The connect component of the identity in the group G{M) of algebraic 
symmetries of a n-dimensional massive Frobenius manifold is a n-dimensional commuta- 
tive Lie group that acts localy transitively on M . 

This is a reformulation of the main lemma of [39]. 

Action of the group of algebraic symmetries provides a new affine structure on a 
massive Frobenius manifold. The structure tensor cji- is constant in this affine structure. 
From the main lemma the following statement follows. 



f^{u ■ v) = f^{u) ■ f^{v) 



(1.13) 



[w, u ■ v] — [w.,u\ ■ V -\- [w, v] ■ u 



(1.14) 
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Theorem 1.1. [39] On a massive Probenius manifold local coordinates u^,...u^ exist 
such that the multiplication law of vector fields in these coordinates has the form 

di ■ dj = 5ijdi, (1-15) 

where di = d/du^. The invariant metric r] in these coordinates has a diagonal form 

n 

r^apdrdtl^ = Y,m{u){du'f (1.16) 



satisfying the equations 



diJ2vii{u)du') = 0, (1.17a) 



i=l 
n 



Conversely, for aflat diagonal metric with the properties (1-17) andt°' = t°'{u), a = 1, ...,n 
being the fiat coordinates for the metric the formulae 

i=l 
i=l 

determine (localy) a massive Frobenius manifold. 

The above coordinates u^, ... , m"^ on a massive Frobenius manifold are determined 
uniquely up to permutations and shifts. They are called canonical coordinates on the 
massive Frobenius manifold M. The tensor c of structure constants in these coordinates 
has the following canonical constant form 

4 = %<^^. (1.19) 

The canonical coordinates can be found by solving an overdetermined system of differ- 
ential equations 

For massive conformal invariant Frobenius manifolds (see the next section) they can be 
found in a pure algebraic way (the Proposition 2.4 below). 
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To complete local classification of massive TFT one has to classify flat diagonal metrics 
with the properties (1.17). This class of metrics was studied by Darboux [33] and Egoroff. 
Following Darboux, I will call them Egoroff metrics. Vanishing of the curvature of these 
metrics can be written in the form of the following system of PDE [Darboux - Egoroff 
system) for the rotation coefficients 







(1.20) 


dklij = liklkj, h j, k are distinct, 




(1.21a) 


n 
k=l 




(1.216) 


Iji = lij- 




(1.21c) 


It is interesting that the same equations (for even n) arise in the calculation [34] of multi- 
point correlators in impenetrable Bose-gas, see Appendix to [45] . 

Integrability of the Darboux - Egoroff system was observed in [23]. It essentialy 
coincides with the "pure imaginary reduction" of the n-wave system (see [26, 25]). This 
can be represented as the compatibility conditions of the following linear system (depending 
on a spectral parameter A) 

djij^i = ^ijipj, i^j (1.22a) 


n 




(1.226) 


To complete local classification of massive Frobenius manifolds one first should apply 
an appropriate version of inverse spectral transform (1ST) to solve the Darboux - Egoroff 
system (1.21). Below I will give an example of 1ST for the important case of self-similar 
solutions of (1.21) (so called topological conformal field theories [8, 51]). To find the metric 
(1.16) and flat coordinates t"' — t^{u) for a given solution iij{u) one has to fix a basis 
ipiaiu), a = 1, ...n in the space of solutions of the system (1.22) for A = 






(1.23a) 


fc 




(1.236) 


a = 1, n. Then we put 

r]ii{u) = ipai'^), 




(1.24a) 


n 

Va/3 = Yi^ia{u)i^if3{u), 
i=l 




(1.246) 


dta 

= V'a(w)V'ia(w), 




(1.24c) 
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n 



{1.2M) 



These formulae complete local classification of complex massive Frobenius manifolds. They 
are parametrized (localy) by n arbitrary functions of 1 variable (the parametrization of 
solutions of the Darboux - Egoroff system) and also by n complex parameters because of 
the ambiguity in the choice of solutions ifjn in the formulae (1.24). 

To classify real Frobenius manifolds one should apply 1ST to various real forms of the 
Darboux - Egoroff system. We will not do it here (see [27] for discussion of real forms of 
the system (1.21) in algebraic-geometry 1ST). 

Global topology of massive Frobenius manifolds is rather poor. We say that a n- 
dimensional manifold M admits S'^j.-structure if the structure group of the tangent bundle 
TM can be reduced to the symmetric group S^- An atlas of coordinates charts on M is 
compatible with the given 5'n-structure if differentials of the transition functions are the 
correspondent elements of Sn (in the standard n-dimensional representation). Globaly 
a 5'^-manifold M with a compatible atlas is determined by an affine representation of 
TTi{M) —>■ Sn ^ Affnj i-G. the transition functions have the form 



for cr, a' G Sn- As an example of S'„-manifold one can have in mind the space of all 
polynomials M = {P{u) = u"' + aiu"'~^ + ... + ai, an G C} without multiple roots. 
Compatible coordinates are the roots of P{u). The transition functions (1.25) are given 
by the standard n-dimensional representation of the braid group 7ri(M) = Bn- 

The Darboirx - Egoroff system is well-defined on any jS^-manifold M with a marked 
compatible atlas. To obtain a massive Frobenius structure on M one should find a solution 
^ij{u) being covariant with respect to transformations of the form (1.25). This "boundary 
value problem" seems to be more complicated. 

In all the examples (below) of massive TFT the coupling space M (massive Frobenius 
manifold) can be extended by adding certain locus Mging (at least of real codimension 2). 
The structure of Frobenius manifold can be extended on M = M U Mging but the algebra 
structure on the tangent spaces T^M for x G Mging has nilpotcnts. The fiat metric rja/s is 
extended on M without degeneration. So M is still a localy Euclidean manifold. 

Remark. The notion of Frobenius manifolfd admits algebraic formalization in terms 
of the ring of functions on a manifold. More precisely, let -R be a commutative associative 
algebra with a unity over a field k of characteristics ^ 2. We are interesting in structures 
of Frobenius algebra over R in the i?-module of A;-derivations Der{R) (i.e. u{k) = for 
K & k,u & Der{R)) satisfying 



(1.25a) 




(1.256) 



V«(A)V„(A) - V„(A)V^(A) = V[^,„](A) identicaly in A 
for Vu(A)i; = Vui) + Xu ■ v, 
V„e = for all u G Der{R) 



(1.26a) 



(1.266) 
(1.26c) 
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where e is the unity of the Frobenius algebra Der{R). Non-degenerateness of the symmetric 
inner product 

< , >: Der{R) x Der{R) R 

means that it provides an isomorphism Y{om.ii{Der{R)^R) Der{R). I recall that the 
covariant derivative is a derivation VuV G Der{R) defined for any u, v & Der{R) being 
determined from the equation 

< V^v, w >= 

^[u < v,w > +v < w,u > —w < u,v > + < [u,v],w > + < [w,u],v > + < [w,v],u>] 

(1.27) 

for any w G Der{R) (here [ , ] denotes the commutator of derivations). Note that the 
notion of infinitesimal algebraic symmetry also can be algebraicaly formalized in a similar 
way. It would be interesting to find a pure algebraic version of the theorem 1.1 . This 
could give an algebraic approach to the problem of classification of Frobenius manifolds. 

We consider in conclusion of this section a closure of the class of massive Frobenius 
manifolds as the set of all Frobenius manifolds with n-dimensional commutative group 
of algebraic symmetries. Let ^ be a fixed n-dimensional Frobenius algebra with struc- 
ture constants c^j and an invariant inner nondegenerate inner product e = (e^). Let us 
introduce matrices 

Ci = (4). (1.28) 
An analogue of the Darboux - Egoroff system (1.21) for an operator- valued function 

^iu):A^A, 7=(7^(t,)), u={u\...,u^) (1.29a) 

(an analogue of the rotation coefficients) where the operator 7 is symmetric with respect 
to e, 

€7 = 7^6 (1.296) 

has the form 

[Ci, dj^] - [Cj, da] + [[Ci, 7] , [Cj, 7]] = 0, z, j = 1, . . . , n, (1.30) 

di = d/du^. This is an integrable system with the Lax representation 

ai* = *(AQ + [Q,7]), z=l,...,n. (1.31) 

It is convenient to consider ^ = {il^iiu), . . . , '0„(u)) as a function with values in the dual 
space A* . Note that A* also is a Frobenius algebra with the structure constants = c^^e*-' 
and the invariant inner product < , >* determined by (e*-^) = {eij)~^ . 
Let '^oi{u), a = l, ... ,n be a basis of solutions of (1.31) for A = 

a*« = *«[Ci,7], a=l,...,n (1.32a) 

such that the vector \E'i(u) is invertible in A* . We put 

Vap =< *a(w), */?(«) >* (1.326) 
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grad„t« = (1.32c) 

CM) = (i.32d) 

Theorem 1.2. Formulae (1.32) for arbitrary Probenius algebra A localy parametrize 
all Frobenius manifolds with n-dimensional commutative group of algebraic symmetries. 

Considering m as a vector in A and = \&i • as a linear function on A one obtains 
the following analogue of Egoroff metrics (on A) 

ds'^ = ^l{du-du). (1.33) 

2. Conformal invariant Frobenius manifolds and isomonodromy deforma- 
tions. 

Definition 2.1. A diffeomorphism / : M — > M is called conformal symmetry if 

f^{u-v)^ tJi)f^{u)- f^(v) (2.1a) 

< f*{u),f4v) >= ii"} <u,v > (2.1b) 

f.{e) = fi}e (2.1c) 

for some functions //^ /Uj, /U^-. A Frobenius manifold M is called conformal invariant if it 
admits a one-parameter group of conformal symmetries /'-'^•' such that the tensors f*^\c), 

f*~\'n): f*^\^) determine on M a Frobenius structure for any r. 

Let V be the generator of the one-parameter group of conformal symmetries on a 
conformal invariant Frobenius manifolds. 

Proposition 2.1. On a massive conformal invariant Frobenius manifold an action 
of the one-parameter group of conformal symmetries is generated by the field 



v = J2'^'9i (2.2) 



i=l 

(modulo obvious transformations v ^ av + be for constant a and b). It acts on the tensors 
c, r], e by the following formulae 

C^c — c (2.3a) 

L^e = -e (2.36) 

jC^r] = (2 - d)r] (2.3c) 

where d is a constant. 

Here Cy means the Lie derivative along the vector field v. 

Corollary. For a massive conformal invariant Frobenius manifold the rotation coef- 
ficients Jij{u) satisfy the similarity condition 

^ij{cu) = c~^7y(tt) (2.4a) 
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or, equivalently 

n 

J2^''^kl^j{u) ^ -'y^J{u). (2.46) 

k=l 

For n — 2 the similarity reduction (2.4) of the Darboiix ~ Egoroff system can be solved 
immediately: 

1 

2 — 

For the first nontrivial case n = 3 the system (1.21), (2.4) reads 

r'l = FaFs (2.6a) 

{zT2y = -F1F3 (2.66) 
((^-l)F3)' = FiF2 (2.6c) 



where 



lijiu) — —5 ^Ffc(^), i, J, k are distinct (2.7a) 



— 



{2.7b) 



It has an obvious first integral 

F^ + {zr2f + {{z - 1)F3)^ = const. (2.8) 

Using this integral one can reduce [30] the system (2.6) to a particular case of the Painleve- 
VI equation. 

For n > 3 the system (1.21), (2.4) can be considered as a high-order analogue of the 
Painlcvc-VI. To find solutions of this system one can use an appropriate version of 1ST: 
the so-called method of isomonodromy deformations [31]. This gives parametrization of 
solutions of the system (1.21), (2.4) by monodromy data of the following system of linear 
ODE with rational coefficients: 

A^ = (AC/-[C/,7])V'. (2.9a) 

Here 

t/ = diag(w\...,ti"), (2.96) 
7=(7ii(«))- (2.9c) 

Solutions of this linear ODE have some monodromy properties, i.e. they are multivalued 
functions in the complex A-planc. 

Proposition 2.2. The monodromy transformations of solutions of the system (2.9) 
do not depend on the parameters u iff the matrix ^ij{u) is a solution of the system (1.21), 

(24). 
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The linear system (2.9) has two singular points: a regular singularity at A = and 
an irregular one at A = oo. Monodromy transformations of solutions of the system near 
A = have the form 

^ ^ exp(-27ri[t/, 7])'!/'. (2.10) 

So the eigenvalues of the matrix [t/, 7] are first integrals of the system (1.21), (2.4). These 
generalise the first integral (2.8). Monodromy at infinity is determined by a n x n Stokes 
matrix S (see [31] for details). The diagonal terms of S equal 1; n{n — l)/2 of other entries 
of the matrix S vanish. Other matrix elements of S can be used as local parameters 
of massive conformal-invariant Frobenius manifolds (just n{n — l)/2 arbitrary complex 
parameters; one should add one more parameter: a norming constant of a solution t/jniu) 
in (1.24) being an eigenvector of the matrix [[/, 7]). The monodromy at A = can be 
expressed via S using cyclic relations (see [39]). If the Stokes matrix S is sufficiently close 
to the unity matrix then the inverse problem of the monodromy theory (i.e., to determine 
the coefficients of the linear operator (2.9) from the given monodromy data) always is 
solvable. The solution can be obtained by solving linear integral equations [39]. 

Let us assume that the monodromy of the operator (2.9) in the origin is semisimple. 
That means that the matrix [U,'y] has pairwise different eigenvalues fii, ... Hn- Let us 
order them in such a way that 

//a + //n-a+l = 0. (2.11) 

Proposition 2.3. Flat coordinates on a massive conformul invariant Frobenius man- 
ifold with semisimple monodromy of (2.9) at X = can be chosen in such a way that the 
generator v of conformal symmetries has the form 

V = J](l - q^)t^d^ (2.12) 

for 

ga = yUl - /Ua (2.13a) 

where Ha are the eigenvalues of the matrix [f/, 7] ordered as in (2.11). 

In other words, the tensors c, ry, e should be conformal covariant with respect to the 
following transformations 

r ^ c^-««t°= (2.14a) 
:7 ,7 (2.146) 



c 



a/3 " ^ ^a/3 

VaP ^ C^^+^^-^rya/? (2.14c) 

d^Qn^ 2//1 (2.136) 

e ^ c~^e. {2.Ud) 

The equation (2.14c) means that rjap 7^ only for Qa + = d. 

The numbers Qce are called charges of the TCFT model, d is called dimension of the 
model. For topological sigma-models it coincides with complex dimension of the target- 
space. Scaling laws (2.14) were obtained in [8] using the assumption that the TCFT model 



where 

is the same as in (2.3c), 
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is obtained by twisting of a N=2 supersymmetric model of QFT. These imply superselection 
rules for trcc-lcvcl correlators in the conformal point t = (the stationary point of the 
field V (2.3)). In our approach the scaling laws follow from simple symmetry assumption 
on the Frobenius manifold. 
Summarizing we obtain 

Theorem 2.1. All massive conformal invariant Frobenius manifolds are parametrized 
by monodromy data of the linear operator 

A = Xdx- XU + M{u) (2.15a) 

C/ = diag(w\...,w"). (2.156) 
= -M. (2.15c) 

Manifolds with semisimple monodromy at the origin A = form a [lillLJl _|_ l'\-parameter 
family. The free energy F{t) of such a Frobenius manifold can be expressed via quadratures 
of a high-order analogue of the Painleve-VI transcendents, i.e. solutions of the equations 
of isomonodromy deformations of (2.15). 

For nonresonant conformal invariant Frobenius manifolds (see (3.19) below) with a 
semisimple monodromy at the origin the structure functions c^^ {t) can be expressed alge- 
braicaly (i.e. without quadratures) via the above high-order analogue of the Painleve-VI 
transcendents. Also one has 

Proposition 2.4. The canonical coordinates , . . . .u^ on a massive conformal in- 
variant Frobenius manifold coincide with eigenvalues of the matrix 

U = (U'^it)) = ((1 + q^- q,)F;{t)) (2.16a) 
Fj{t) = r^'^'dpd,F{t). (2.166) 

It would be interesting to understand a physical sense of the operator U for TCFT 
models. 

Remark. We saw that monodromy is an important invariant of a massive conformal 
invariant Frobenius manifold. It can be defined also for arbitrary conformal invariant 
Frobenius manifold by considering the linear operator 

K=\dx-\U + M, (2, 17a) 

where 

M = {MJ) = {qp5}), (2.176) 

the matrjx U has the form (2.16). WDVV equations determine isomonodromy deforma- 
tions of A. 

Monodromy properties of eigenfunctions of A near irregular singularity A = oo (i.e. 
Stokes matrices) strongly depend on algebraic structure of the multiplication on TM. 
These Stokes matrices are constrained by cyclic relations since monodromy near the origin 
A = is fixed by the given charges qa- An advantage of the isomonodromy problem 
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(2.15) for massive Frobenius manifolds is in universality (independence on the charges; the 
charges can be expressed via an arbitrary Stokes matrix of (2.15)). Note that a basis of 
common eigenfunctions of A 

AC = (2.18a) 

and (1-12) has the form 

^I3{t, X) = d/jhait, X), K = d — qa, for any Q! = 1, . . . , n (2.186) 

where the solutions hce{t, A) of (3.5) are normalized by (3.6). 

3. Coupling to gravity. Systems of hydro dynamic type: their Hamiltonian 
formalism, solutions, and r-functions. 

Let us fix a Frobenius manifold (i.e. a solution of the WDVV equations. Considering 
this as the primary free energy of the matter sector of a 2D TFT model, let us try to calcu- 
late the tree- level (i.e., the zero-genus) approximation of the complete model obtained by 
coupling of the matter sector to topological gravity. The idea to use hierarchies of Hamil- 
tonian systems of hydrodynamic type for such a calculation was proposed by E.Witten 
[46] for the case of topological sigma-models. An advantage of my approach is in effective 
construction of these hierarchies for any solution of WDVV. The tree- level free energy of 
the model will be identified with r-function of a particular solution of the hierarchy. For 
a TCFT-model (i.e. for a conformal invariant Frobenius manifold) the hierarchy carries 
a bihamiltonian structure under a non-resonance assumption for charges and dimension 
of the model (this bihamiltonian structure was constructed in [39] for the case of massive 
perturbations of a TCFT model; here I generalize it for an arbitrary TCFT model). This 
gives an answer to a question of [46] (see p. 283). As it was mentioned in the Introduction, 
the bihamiltonian structure could be useful for calculation higher genus corrections . 

So let c^^(t), r]cei3 be a solution of WDVV, t = (t^, . . . , f^). I will construct a hierarchy 
of the first order PDE systems linear in derivatives {systems of hydrodynamic type) for 
functions t°'{T)^ T is an infinite vector 

T=(T«'^), a = l, ...,n, p = 0, 1, . . . ; T^'^ = X, 

dTo.,A^ = c^o.,p)^^{t)dxt^ (3.1a) 

for some matrices of coefficients C(q, p)^(t). The marked variable X = T^'^ usualy is called 
cosmological constant. 

I will consider the equations (3.1) as dynamical systems (for any (ct,p)) on the space 
of functions t = t{X) with values in the Frobenius manifold M. 

A. Construction of the systems. I define a Poisson bracket on the space of functions 
t — t{X) (i.e. on the loop space C{M)) by the formula 

{t'^{X),t^{Y)} = r]'^^5'{X-Y). (3.2) 

All the systems (3.1a) have hamiltonian form 

dTc.,vt^ = {t^{X),H^,p} (3.16) 
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with the Hamiltonians of the form 

Ha,p = j K,p+i{t{X))dX. (3.3) 
The generating functions of densities of the Hamiltonians 

oo 

/i«(t,A) = J]/i«,p(t)A^ a = l,...,n (3.4) 

p=0 

coincide with the flat coordinates of the perturbed connection V(A) (see (1.10)). That 
means that they are determined by the system (cf. (1.12)) 

df3djha{t,X) = Xc''^^{t)d^ha{t,X). (3.5) 

This gives simple recurrence relations for the densities ha,p- Solutions of (3.5) can be 
normalized in such a way that 

/la(t, 0) = ta = Vapt'^, (3.6a) 

< Vh^{t, X),Vh^{t, -A) >= rj^0. (3.66) 

Here V is the gradient (in t). It can be shown that the Hamiltonians (3.3) are in involution. 
So all the systems of the hierarchy (3.1) commute pairwise. 

B. Specification of a solution t — t{T). The hierarchy (3.1) admits an obvious scaling 
group 

T^,p^cT'''P, t^t. (3.7) 
Let us take the nonconstant invariant solution for the symmetry 

(driA - ^T"'faT«.p)t(T) = (3.8) 

(I identify T^'^ and X. So the variable X is supressed in the formulae.) This solution can 
be found without quadratures from a fixed point equation for the gradient map 

t = V$t(^), (3.9) 
$T(t) = J]T"'P/i«,p(t). (3.10) 

a,p 

It can be proved existence and uniqueness of such a fixed point for sufficiently small T"'^ 
for p > (more precisely, in the domain: T'*'^ are arbitrary, T^'^ = o(l), T'^'P = o{T^'^) 
for p > 0). 

C. T-function. Let us define coefficients V(^cx,p),{P,q)i^) ^^om the expansion 

oo 

(A + //)-!(< Vha{t,X),Vhp{t,ii) > -r]^j3) = ^ia,pUi3,q){t)XPn'' = Vo,p{t,X,n). 

p,q=0 

(3.11) 
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The infinite matrix of coefficients V(a,p),(/3,q)(^) has a simple meaning: it is the energy- 
momentum tensor of the commutative Hamiltonian hierarchy (3.1). That means that a 
matrix entry V(a,p),(/3,q)(^) is the density of flux of the Hamiltonian H^ p along the flow 

dT0,iha,p+l{t) = dxV(^a,p),{l3,q){t). (3-12) 

Then 

(3.13) 

Remark. More general family of solutions of (3.1) has the form 

V[$T(t)-$To(0] = (3.14) 

for arbitrary constant vector Tq = Tq'^. For massive Frobenius manifolds these form a 
dense subset in the space of all solutions of (3.1) (see [22, 23, 39]). Formally they can 
be obtained from the solution (3.9) by a shift of the arguments T"'^. r-function of the 
solution (3.14) can be formaly obtained from (3.13) by the same shift. For the example of 
topological gravity [3, 46] such a shift is just the operation that relates the tree-level free 
energies of the topological phase of 2D gravity and of the matrix model. It should be taken 
in account that the operation of such a time shift in systems of hydrodynamic type is a 
subtle one: it can pass through a point of gradient catastrophe where derivatives become 
inflnite. The correspondent solution of the KdV hierarchy has no gradient catastrophes 
but oscillating zones arise (see [32] for details). 
Theorem 3.1. Let 

^(T) = logT(T), (3.15a) 

< (f>a,p(f>(3,q ■ ■ ■ >0= dT'^,pdTl3,q . . .T{T). (3.156) 

Then the following relations hold 

T{T)\ n =F{t) (3.16a) 

\ J\t°''P=0 for p>0, TO!. 0=40! \ / \ / 

dxT{T) = J2 T"'^aTo,P-iJ^(T) + IrjapT^'^Tf^^^ (3.166) 

< (I>a,p(l>l3,q(l>j,r >0==< 4>cx,p-l4>X,0 >0 V^'' < 4>l^,0(Pl3,q4>'r,r >0 • (3.16c) 

Let me establish now a 1-1 correspondence between the statements of the theorem 
and the standard terminology of QFT. In a complete model of 2D TFT (i.e. a matter 
sector coupled to topological gravity) there are infinite number of operators. They usualy 
are denoted by 4'ce,p oi (Tp{4>ce)- The operators (pce^ can be identified with the primary 
operators (p^', the operators (j)a,p for p > are called gravitational descendants of ^q,. 
Respectively one has infinite number of coupling constants T"'^. The formula (3.15a) 
expresses the tree-level (i.e. genus zero) partition function of the model of 2D TFT via 
logarythm of the r-function (3.13). Equation (3.15b) is the standard relation between the 
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correlators (of genus zero) in the model and the free energy. Equation (3.16a) manifests 
that before coupling to gravity the partition function (3.15a) coincides with the primary 
partition function of the given matter sector. Equation (3.16b) is the string equation for 
the free energy [3, 4, 8, 46]. And equations (3.16c) coincide with the genus zero recursion 
relations for correlators of a TFT [4, 46] . 
Particularly, from (3.15) one obtains 

< (Pa,p(p/3,q >0= V^a,p),{p,q){t{T)), (3.17a) 

< (f>a,p(f>i,o >o= K,p{t{T)), (3.176) 

< (t>a,p<t>P,q4>^,r >0 = < V/la,p • V/l/3,g • V/l^,r, [e - V/la,p-l] > • (3.17c) 

The second factor of the inner product in the r.h.s. of (3.17c) is an invertible element (in 
the Frobenius algebra of vector fields on M) for sufficiently small T"'^, p > 0. Prom the 
last formula one obtains 

Proposition 3.1. The coefficients 

^pZcP (^) = r]^^dTc^,vdT,,.dT„v logT(T) (3.18) 

for any p and any T are structure constants of a commutative associative algebra with the 
invariant inner product rjafs- 

As a rule such an algebra has no unity. 

In fact the Proposition holds also for a r-function of an arbitrary solution of the form 
(3.14). 

We see that the hierarchy (3.1) determines a family of Backlund transforms of the 
WDVV equation (0.1) 

F{t) ^ Fit), 
F = logT, i'^ = T''^P 

for a fixed p and for arbitrary r-function of (3.1). So it is natural to consider equations of 
the hierarchy as Lie - Backlund symmetries of WDW. 

Up to now I even did not use the scaling invariance (2.14). It turns out that this gives 
rise to a bihamiltonian structure of the hierarchy (3.1). 

Let us consider a conformal invariant Frobenius manifold, i.e. a TCFT model with 
charges Qa and dimension d. We say that a pair a,p is resonant if 

^1-q^+p^O. (3.19) 

Here p is a nonnegative integer. The TCFT model is nonresonant if all pairs a,p are 
nonresonant. For example, models satisfying the inequalities 

= Ql < < . . . < gn = < 1 (3.20) 

all are nonresonant. 



23 



Theorem 3.2. 1) For a conformal invariant Frobenius manifold with charges Qq, and 
dimension d the formula 



{t-{X), = [(^ - qa)F-^t{X)) + - qf,)F-^{t{Y))]6\X - Y) (3.21) 

F"^(t) = jf^'jf'^'dc,>dp>F{t) 



determines a Poisson bracket compatible with the Poisson bracket (3.2). 2) For a nonres- 
onant TCFT model all the equations of the hierarchy (3.1) are Hamiltonian equations also 
with respect to the Poisson bracket (3.21). 

The nonresonancy condition is essential: equations (3.1) with resonant numbers {a,p) 
do not admit another Poisson structure. 

Remark. According to the theory [18-21] of Poisson brackets of hydrodynamic type 
any such a bracket is determined by a flat Riemannian (or pseudo-Riemannian) metric 
da/sit) on the target space M (more precisely, one needs a metric g°'^{t) on the cotangent 
bundle to M). In our case the target space is the Frobenius manifold M. The first 
Poisson structure (3.2) is determined by the metric being specified by the double-point 
correlators rja/B- The second fiat metric for the Poisson bracket (3.21) on a conformal 
invariant Frobenius manifold M has the following geometrical interpretation. Let ui and 
u>2 be two 2-forms on M. We can multiply them uji, uj2^ •<^2 using the multiplication 
of tangent vectors and the isomorphism r] between tangent and cotangent spaces. Then 
the new inner product < , >i is defined by the formula 



Here i^ is the operator of contraction with the vector field v (the generator of confor- 
mal symmetries (2.3)). The metric (3.22) can be degenerate. The theorem states that, 
nevertheless, the Jacobi identity for the Poisson bracket (3.21) holds. 

Main examples of solutions of WDVV and of corresponding hierarchies will be given 
in the next section. Here I will consider the simplest class of examples where c^^ does not 
depend on t. They form structure constants of a Frobenius algebra A with an invariant 
inner product < , > {r]ai3 in a basis ei = 1, . . . , e^)- Let 



<uji,uj2 >i= i„(a;i • UJ2). 



(3.22) 



t = t"ea e A. 



(3.23) 



The linear system (3.5) can be solved easily: 

ha{t, A) = < Ca, e 



At 



- 1 > . 



This gives the following form of the hierarchy (3.1) 

p\ 



(3.24) 



The solution (3.9) is specified as the fixed point 



G{t) = t. 



(3.25a) 
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Here I introduce A-valued coupling constants 

Tp = T"'fe« e A P=0, 1,.... (3.26) 
The solution of (3.25) has the well-known form 

t = G(G(G(...))) (3.27) 
(infinite number of iterations). The r- function of the solution (3.27) has the form 



^ 6^ ' ^ ^ (p+2)p! ^2^ + 5 + 



(3.28) 



For the tree-level correlation functions of a TFT-model with constant primary correlators 
one immediately obtains 

^ ^ ^ ^ <ec.ef3,tP+i+^ > 

< (p + g+l)p!g! ' ('-''"^ 

< <t>a,p<t>P,q<i>^,r >0= -T-j-J < eaC/^e^, — — T^r^ > . (3.296) 

2^s>l (s-1)! 

For n = 1 the formulae (3.29) give the tree-level correlators of the topological gravity (see 
[3, 46]). For n = 24 one obtains the tree-level correlators of the topological sigma- model 
with a K3-surface as the target space. Here the algebra A = H* {K3) is a graded one: it 
has a basis P, Qi, Q22, R of degrees 0, 1 (all the Q's) and 2 resp. The multiplication 
has the form 

P is the unity, QiQj = r]ijR, QiR ^ ^ (3.30) 

for a nondegenerate symmetric matrix rjij. The scalar product (the intersection number) 
has the form 

Vpr = 1, VQiQj = Vij- 

Let us consider now the second hamiltonian structure (3.21). I start with the most 
elementary case n — 1 (the pure gravity). Let me redenote the coupling constant 

u = t^. 

The Poisson bracket (3.21) for this case reads 

MX), u{Y)h = l{u{X) + u{Y))5'{X - Y). (3.31) 
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This is nothing but the Lie - Poisson bracket on the dual space to the Lie algebra of 
one-dimensional vector fields. 

For arbitrary graded Frobenius algebra A the Poisson bracket (3.21) also is linear in 
the coordinates 

{t-{X),t''{Y)}, = [(^ - qa)cft^X) + - q^)cft'r{Y)]6'{X - Y). (3.32) 

It determines therefore a structure of an infinite dimensional Lie algebra on the loop 
space C{A*) where A* is the dual space to the graded Frobenius algebra A. Theory of 
linear Poisson brackets of hydrodynamic type and of corresponding infinite dimensional 
Lie algebras was constructed in [34] (see also [18]). But the class of examples (3.32) is a 
new one. Note that the case A = H* {K2>) is a nonresonant one. 

Let us come back to the general (i.e. nonlinear) case of a TCFT model. I will assume 
that the charges and the dimension are ordered in such a way that 

= < < • • • < qn-i <qn^d. (3.33) 
Then from (3.21) one obtains 

{r(x),r(r)}i = ^^{e{x) + e{Y))5\x-Y). (3.34) 

Since 

{t-{X\ r(F)}i = [(^ - q^)r{X) + l^t"(F)]5'(X - F), (3.35) 
the functional 

P = J t''{X)dX (3.36) 

generates spatial translations. We see that for d ^ 1 the Poisson bracket (3.21) can be 
considered as a nonlinear extension of the Lie algebra of one-dimensional vector fields. 
An interesting question is to find an analogue of the Gelfand - Fuchs cocycle for this 
bracket. I found such a cocycle for a more particular class of TCFT models. Wc say that 
a TCFT- model is graded if for any t the Frobenius algebra c^^(t), r}ai3 is graded. 
Theorem 3.3. For a graded TCFT-model the formula 

{r{X),t^{Y)}^ = {t''{X),t^{Y)}i + eWW^S"'iX - Y) (3.37) 

determines a Poisson bracket compatible with (3.2) and (3.21) for arbitrary (the central 
charge). For a generic graded TCFT model this is the only one deformation of the Poisson 
bracket (3.21) proportional to 5"'{X - Y). 

For n=l (3.37) determines nothing but the Lie - Poisson bracket on the dual space 
to the Virasoro algebra 

{u{X), u{Y)}^ = ^[u{X) + u{Y)]5'{X -Y) + e''5"\X - Y) (3.38) 
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(the second Poisson structure of the KdV hierarchy). For n > 1 and constant primary- 
correlators (i.e. for a constant graded Frobenius algebra A) the Poisson bracket (3.37) can 
be considered as a vector-valued extension (for d ^ 1) of the Virasoro. 

Graded TCFT models occur as the topological sigma-models with a Calabi - Yau 
manifold of (complex) dimension d as the target space [2, 46, 57]. They are nonresonant 
for even d. Particularly, for d = 2 one obtains the K3-models where the primary correlators 
are constant. For d > 2 these are not constant because of instanton corrections [46, 47, 
57]. As it was explained in [57], finding of these primary correlators for the Calabi - Yau 
models (and, therefore, graded solutions of WDVV) could be a crucial point in solving the 
problem of mirror symmetry. 

The compatible pair of the Poisson brackets (3.2) and (3.37) generates an intcgrable 
hierarchy of PDE for a non-resonant graded TCFT using the standard machinery of the 
bihamiltonian formalism [52] 

= {t^(X),^«,p} = {tf'{X),H^,p_^}\. (3.39) 
Here the Hamiltonians have the form 

Ha,p = J ha,p+idX, (3.40a) 

ha,p+i = [^^-^-qa+p]~^ha,pMt) + ^'^^ha,p+i{i^dxt,...,d^t-,e^) (3.406) 

where Aha,p+i are some polynomials determined by (3.39). They are graded-homogeneous 
of degree 2 where degd^t — k, dege = — 1. The small dispersion parameter e also plays 
the role of the string coupling constant. It is clear that the hierarchy (3.1) is the zero- 
dispersion limit of this hierarchy. For n = 1 using the pair (3.2) and (3.38) one immediately 
obtains the KdV hierarchy. Note that this describes the topological gravity. It would be 
interesting to investigate relation of the hierarchies determined by the pair (3.2) and (3.37) 
to a nonperturbative (i.e., for all genera) description of the Calabi - Yau models (especially, 
of the K3 models) coupled to gravity. For a model with constant correlators (for a graded 
Frobenius algebra A) the first nontrivial equations of the hierarchy are 

^T".it = e^ttx + 3 _ ^ ^a^ntxxx- (3-41) 

For non-graded TCFT models it could be of interest to find nonlinear analogues of 
the cocycle (3.37). These should be differential geometric Poisson brackets of the third 
order [58, 18] of the form 

{r{x),t\Y)}, = {r{x),t^{Y)}^+ 

e^{g''P{t{X))5"\X - F) + hf {t{X))tl5" {X - Y) + 
[ff{t{X))tl^ + hf,{t{X))tlt'^]5'{X - y)+ 
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\pf{t)tlxx + <f + rf,^{t)tlt'xt),nx - Y)}. (3.42) 

I recall (see [58, 18]) that the form (3.42) of the Poisson bracket should be invariant 
with respect to nonlinear changes of coordinates in the manifold M. This implies that the 
leading term g'^^{t) transforms like a metric (may be, degenerate) on the cotangent bundle 
T^M^ b'^^{t) are contravariant components of a connection on M etc. The Poisson bracket 
(3.42) is assumed to be compatible with (3.2). Then the compatible pair (3.2), (3.42) of 
the Poisson brackets generates an integrable hierarchy of the same structure (3.39), (3.40). 

4. Examples. 

I start with the most elementary examples of solutions of WDVV for n = 2. Only 
massive solutions are of interest here (a 2-dimensional nilpotent Frobenius algebra has 
no nontrivial deformations). The Darboux - Egoroff equations in this case are linear. I 
consider only TCFT case (the similarity reduction of WDVV) . Let us redenote the coupling 
constants 

= u, ^ p. (4.1) 
For 1 the primary free energy F has the form 

« = ^ (4.3) 

g is an arbitrary constant. The second term in the formula for the free energy should be 
understood as 

The linear system (3.5) can be solved via Bessel functions [39]. Let me give an example of 
equations of the hierarchy (3.1) (the T = T^'^-flow) 

ut + uux + gp^'px = (4.4a) 

Pt + ipu)x = 0. (4.46) 

These are the equations of isentropic motion of onc-dimensional fluid with the dependence 
of the pressure on the density of the form p = -^^p"^^- The Poisson structure (3.2) for 
these equations was proposed in [37]. For a = (equivalently d = —1) the system coincides 
with the equations of waves on shallow water (the dispersionless limit [59] of the nonlinear 
Schrodinger equation (NLS)). 

For d = 1 the primary free energy has the form 

F = ]-pu^ + geP. (4.5) 

This coincides with the free energy of the topological sigma-model with CP^ as the target 
space. Note that this can be obtained from the same solution of the Darboux - Egoroff 
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system as the semiclassical limit of the NLS (the case d = —1 above) for different choices 
of the eigenfunction ipu (in the notations of (1.24)). The corresponding T = T^'°-system 
of the hierarchy (3.1) reads 

ut = g{e^)x 
Pt = ux- 

Ehminating u one obtains the long wave limit 

Ptt = g{e^)xx (4.6) 

of the Toda system 

p^^^ = e^"+i - 2e^" + e^"-i . (4.7) 

(The 2-dimensional version of (4.6) was obtained in the formalism of Whitham-type equa- 
tions in [44].) It would be interesting to prove that the nonperturbative free energy of the 
CP^-model coincides with the r-function of the Toda hierarchy. 

Example 2. Topological minimal models. I consider here the A„-series models only. 
The Frobenius manifold M here is the set of all polinomials (Landau - Ginshurg superpo- 
tentials) of the form 

M = {w{p) = p'^+i + aij9"-i + . . . + ani ai, . . . , an e C}. (4.8) 

For any w e M the Frobenius algebra A = Ayj is the algebra of truncated polynomials 

= C\p]/{w'{p) = 0) (4.9) 

(the prime means derivative with respect to p) with the invariant inner product 

</,S>=res,_i<^^. (4.10) 

The algebra A^ is semisimple if the polynomial w'{p) has simple roots. The canonical 
coordinates (1.15) tt^, . . . ,tt" are the critical values of the polynomial w{p) 

It* = w{pi), where w'{pi) — 0, i = 1, . . . ,n. (4-11) 

Let us take the following diagonal metric on M 

n 

Y,r]ii{u){du')\ r^ii{u) = [w"{p,)]-\ (4.12) 
1=1 

It can be proved that this is a flat Egoroff metric on M. The correspondent flat coordinates 
on M have the form 

t"' = j^ resp^oow {p)dp, Q; = l,...,n. (4-13) 
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The metric (4.12) in these coordinates has the constant form 

n 

J2vii{u){du'f = riapdf'dt'^, rjap = Sn+l,a+l3- (4.14) 

i=l 

The ortonormal basis in with respect to this metric consists of the polynomials (/)i(p),... 
, 0n(p) of degrees 0, 1, ... , n — 1 resp. where 

d a 

0a(p) = + , Q! = l,...,n, (4.15) 

Here [ ]_(_ means the polynomial part of the power series in p. This is a TCFT model with 
the charges and dimension 

Oi-1 , n-1 /, 

Qa = — -rr, d = qn = — — • (4.16) 
n + 1 n + 1 

In fact one obtains a n-parameter family of TFT models with the same canonical 
coordinates tt* of the form (4.11) where 

riiiiu) ^ r]ii{u, c) = K'(Pi)]~MX] Ca(t>a{Pi)T, (4.17a) 

Th ~|~ J- 

depending on arbitrary parameters ci, c„. This reflects the ambiguity in the choice of 
the solution ipn in the formulae (1.24). These models are conformal invariant if only one 
of the coefficients is nonzero. 

The corresponding hierarchy of the systems of hydro dynamic type (3.1) coincides with 
the dispersionless limit of the Gelfand - Dickey hierarchy for the scalar Lax operator of 
order n+1. This essentialy follows from [8, 11]. I recall that the Gelfand - Dickey hierarchy 
for an operator 

L = + ai(a;)a"-^ + . . . + an{x) 
d = d/dx 

has the form 

at«,.L = Ca,p[L,[L*+^]+], a=l,...,n, p=0,l,... (4.18) 

for some constants Ca,p- Here [ ]+ denotes differential part of the pseudodifferential oper- 
ator. The dispersionless limit of the hierarchy is defined as follows: one should substitute 

x^ex^X, ^ e^'P = T'^'P (4.19) 

and tend e to zero. The dispersionless limit of r- function of the hierarchy is defined [11, 
53-54, 61] as 

logTdispersionless(7') = Am e"^ logT(et). (4.20) 

e— ►0 
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Modified minimal model (4.17) is related to the same Gelfand - Dickey hierarchy with 
the following modification of the L-operator 

L^L^^c^[L^] + . (4.21) 

The linear equation (3.5) for the minimal model can be solved in the form [39] 

Ti ~\~ la a 

ha{t;X) = vesp^oow^ iFi{l;l + Xw{p))dp. (4.22) 

a n + 1 

Here iFi(a; c; z) is the Kummer (or confluent hypergeometric) function [35] 

^^M=J:^^'^, (4.23a) 

{a)m = a{a + l) ...{a + m-1). (4.236) 
The generating function (3.11) has the form 

Vai3{t;X,n) = {X + ii)-^[r]''''{TeSp=^w^~\Fi{l-, Xw{p))(f)f,{p)dp)x 

71/ -L 

{TeSp=ooW^~^iFi(l; iJiw{p))(j)„{p)dp) - r]af3]- (4.24) 

From this one obtains formulae for the r-function. 

Example 3. Mg.^Q^...^„^-models [13, 14]. Let M = Mg.no,...,nm be a moduli space of 
dimension 

n=2g + no + ... + nm + 2m (4.25) 

of sets 

(C; ooo, . . . , oom; w;ko,..., km; ai, . . . ,ag,bi, . . . ,bg) e Mg.no,...,nm (4.26) 

where C is a Riemann surface with marked points ooq, oo^, and a marked meromorphic 
function 

w.C^CP'^, w"^(oo) = oooU...,UoOm (4.27) 

having a degree + 1 near the point oo^, and a marked symplectic basis ai,...,ag, 
bij . . . ,bg G Hi{C, Z), and marked branches of roots of w near ooq, oo^ of the orders 
no + 1, Um + l resp., 

fcJ»^+i(P) = w{P), P near oo^. (4.28) 
(This is a connect manifold as it follows from [56].) We need the critical values of w 

= w{Pj), dw\p. = 0, j = 1, . . . , n (4.29) 

(i.e. the ramification points of the Riemann surface (4.27)) to be local coordinates in open 
domains in M where 

u' ^ for i 7^ j (4.30) 
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(due to the Riemann existence theorem). Another assumption is that the one-dimensional 
affine group acts on M as 

(C; ooo, • • • , oo^; w; . . .) ^ (C; ooq, • • • , oom] aw + b; . . .) (4.31a) 

t-^ au^ + b, i=l,...,n. (4.316) 

Let dp be the normahzed Abehan differential of the second kind on C with a double pole 
at ooo 

dp = dko + regular terms (4.32a) 

dp^O, i^l,...,g. (4.326) 

Using as the canonical coordinates (1.15) I define a flat Egoroff metric on M by the 
formula 

n 

ds'^ = ^r]ii{u){du'f, (4.33a) 

r7n(«) = resp,^. (4.336) 
It can be extended globaly on M. The corresponding flat coordinates are 

f'" = ^ ^ reSoofer~"^^^P: i = 0,...,m, a=l, ...,ni; (4.34a) 

Ui — a + 1 * 

= v.p. / dp = lim ( / dp + ko{Q)), z = 1, . . . , m; (4.346) 

Jooo Q-^ooo Jq 

= —leSoo^wdp, i = 1, . . . , m; (4.34c) 

wdp, i = 1, . . . ,g. (4.34ci) 



f dp, 



The metric (4.33) in the coordinates has the following form 

rii.i;ai;i;0 = — 7 5y 5a+/3,ni +1 (4.35a) 

rii + l 

Vpiqi = 5ij (4.356) 

■qriai = Sij, (4.35c) 

other components of the 77 vanish. The unity vector field is a unit vector along the coor- 
dinate t^'^. 

Proposition 4.1. The flat metric (4-. 35) is well-defined globaly on M and the flat 
coordinates (4-34) are globaly independent analytic functions on M . 

As a consequence we obtain that the moduli space M is an unramified covering over 
a domain in (see [13, 14]). 
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Let us introduce primary differentials on C (or on a universal covering (7 of C \ ooq U 
. . . U oo^) of the form 

= dtA {pdw)^=const (4.36) 
where ^ 

p{P) = [ dp, (4.37a) 

Qo e C, w{Qo) = 0, (4.376) 

is one of the fiat coordinates (4.34). Note that the definition (4.36) of the primary 
differentials can be rewritten as 

(|)^A = -dfA {wdp)p=const (4.38) 

where the multivalued coordinate p on C is defined in (4.37). So w{p) plays the role of the 
Landau - Ginsburg superpotential for the M^.^^ ... „^^-models. More explicitly, (f)ii;a is a 
normalized Abelian differential of the second kind with a pole in ooj, 

(f)^i;a — dk" + regular terms near oOj, 

a 



0; (4.39a) 



(j)pi is a normalized Abelian differential of the second kind on C with a pole only at oOj 
with the principal part of the form 

(j)pi = dw + regular terms near oOj, 

(j)pi = 0; (4.396) 



/ 



(j)qi is a normalized Abelian differential of the third kind with simple poles at ooq and ooj 
with residues -1 and +1 resp.; 

is a normalized multivalued differential on C with increments along the cycles hi of the 
form 

(f)j.i{P + bj) — (f)ri{P) = —Sijdw, 



0; (4.39c) 
(pgi are the basic holomorphic differentials* on C normalized by the condition 

2m5ij. (4.39d) 



* 1-form pdw was used by Novikov and Veselov in their theory of algebro-geomctric 
Poisson brackets [60]. The coordinates are the algebro-geometric action variables of 
[60] . In [60] it also was an important point that derivatives dgi {pdw) are the normalized 
holomorphic differentials. 
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The inner product (4.33) in terms of the primary differentials reads 

,,, = f,e.P.^. (4.40) 

The structure functions CABcif) can be calculated as 

Extension of the Frobenius structure on all the moduli space M is given by the condition 
that the differential 

- c%(l)tcdp 



dw 



(4.42) 



is holomorphic for \w\ < oo. The Frobenius algebra on T^M will be nilpotent for Riemann 
surfaces w : C ^ CP^ with more than double branch points. This is a conformal invariant 
Fr5benius manifold with the dimension 

d = ^ (4.43a) 
no + 1 ^ ^ 

and charges 

?t.« = - (4.436) 

rii + l no + 1 

Qri = Qpi = "° (4.43c) 
no + 1 

Qsi = Qqi = TT- (4.43(i) 

no + 1 

For the particular case g = m = Owe obtain the Frobenius manifolds of the minimal models 
(the previous example). For = 0, m > we obtain models with rational functions as 
superpotentials. 

The generating functions h^A (t; A) (3.4) have the form 

ht,;. {t; A) = -^^reSp=oo,/i;»"iFi(l; 1 + Xw{p))dp. (4.44a) 

a nj + 1 

= v.p. / e^'^dp (4.446) 

J OOn 



'OOo 



hqi = reSooi ^ dp (4.44c) 



hri = <h e^'^dp (4.44cZ) 
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h^. = ^i pe^'^dw. (4.44e) 

Remark. Integrals of the form (4.44) seem to be interesting functions on the moduh 
space of the form Mg.no,---,nm- ^ simplest example of such an integral for a family of elliptic 
curves reads 



/ 

^0 



e^^^^Uz (4.45) 



where p{z) is the Weierstrass function with periods 2(jj, 2uj' . For real negative A a degen- 
eration of the elliptic curve (a; — > oo) reduces (4.45) to the standard probability integral 
Jq°° e^^ dx. So the integral (4.45) is an analogue of the probability integral as a function 
on A and on moduli of the elliptic curve. I recall that dependence on these parameters is 
specified by the equations (3.5), (2.17). 

Gradients of this functions on the moduli space M have the form 



a—rii — l r ^ 
i 1 



dtAhpi = rjtApi - Av.p. 
dfAhqi = resoo^e 



fli -r J- 


(4.46a) 


/ e^-cP.A 


(4.466) 


J OOo 






(4.46c) 


I e^^(l>tA 


(4.46d) 



9t^hsi = ^.j e^^<l>tA. (4.46e) 

The generating function Va/3{t; A, //) of coefficients of the r-function (3.13) can be calculated 
via inner products (w.r.t. the matrix (4.35)) of (4.46). Particularly, for a part of the Hessian 
of the primary free energy F{t) (a function on M) one obtains [13, 14] 

d^F f 

-Tij = - (p 4),i. (4.47) 



This is nothing but the matrix of periods of holomorphic differentials on M. Other second 
derivatives of F also turn out to be certain periods of some Abelian differentials on C. 

Conclusion. WDVV is a universal system of integrable differential equations for 
periods of Abelian differentials on Riemann surfaces. 

I recall that this system is a high-order analogue of the Painleve-VI equation (i.e. 
equations of isomonodromy deformations of (2.15)). To specify the solution of WDVV one 
needs to find the monodromy matrix of the linear operator (2.15) for the eigenfunctions of 
the form (4.44). I will do it in a forthcoming publication. 

We obtain the following picture of "Painleve uniformisation" of the moduli spaces 
Mg-no,...,nm- (1) ^ global system of analytic coordinates on Mg.no,...,nm'i (2) periods of 



35 



Abelian differentials on curves C E Mg.no,---,nm certain high-order Painleve transcen- 
dents as functions of these coordinates. 

Remark. For any Hamiltonian Ha,p of the form (3.3), (4.44) one can construct a 
differential ^A,p on C or on the covering C with singulariries only at the marked infinite 
points such that 

^V,p = resp,^-^, i = l,...,n. (4.48) 

See[13] for an explicit form of these differentials (for m = also see [14]). Using these 
differentials the hierarchy (3.1) can be written in the Flaschka - Forest -McLaughlin form 
[16] 

dTA,pdp = dx^A,p (4.49) 

(derivatives of the differentials are to be calculated with w =const.). 

The matrix V(^A,p),{B,q){t) determines a pairing of these differentials with values in 
functions on the moduli space 

{^A,p, ^B,ci) = V^A,pUB,g) (t) (4.50) 

Particularly, the primary free energy F as a function on M can be written in the form [13, 
14] 

F = --{pdw,pdw). (4-51) 
Note that the differential pdw can be written in the form 

- E + E (4-52) 

where Q^i^^'' is the Abelian differential of the second kind with a pole at ooj of the form 

= ciA;f*+^ + regular terms near ooj. (4.53) 
For the pairing (4.50) one can obtain from [44] the following formula 

ihdw, f2dw) = ^J J^idfidf2 + dhdf2) (4.54) 

where the differentials d and d along the Riemann surface should be understood in the 
distribution sense. The meromorphic differentials fidw and f2d'w on the covering C should 
be considered as piecewise meromorphic differentials on C with jumps on some cuts. 

The corresponding hierarchy (3.1) is obtained by averaging along invariant tori of 
a family of (/-gap solutions of a KdV-type hierarchy related to a matrix operator L of 
the matrix order m + 1 and of orders tiq, rim in d/dx. The example m = (the 
averaged Gelfand - Dickey hierarchy) was considered in more details in [14]. The Poisson 
bracket (3.2) is a result of semiclassical limit (or averaging) [18-21] of the first haniltonian 
structure of the Gelfand - Dickey hierarchy; averaging of the second hamiltonian structure 
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(the classical 1^-algebra) gives the Poisson structure (3.21). Therefore, (3.21) can be 
considered as semiclassical limit of classical VF-algcbras. The corresponding flat metric 
(3.22) on the moduli space Mg.no,...,nm well-defined on a subset of Riemann surfaces 
having w — a. non-ramifying point. 

Also for ^ -h m > one needs to extend the KdV-type hierarchy to obtain (3.1) (see 
[13-14]). To explain the nature of such an extension let us consider the simplest example 
of m = 0, ?io = 1. The moduli space M consists of hyperelliptic curves of genus g with 
marked homology basis 

25+1 

This parametrizes the family of ^-gap solutions of the KdV. The L operator has the well- 
known form 

L = -dl + u. (4.56) 

In real smooth periodic case u{x -\-T) = u{x) the quasimomentum p{w) is defined by the 
formula 

^{x + T,w) = e'^^'^^^^{x,w) (4.57) 
for a solution iIj{x,w) of the equation 

Li/j = wi/j (4.58) 

(the Bloch - Floquet eigenfunction). The differential dp can be extended onto the family 
of all (i.e. quasiperiodic complex meromorphic) f^-gap operators (4.55) as a normalized 
Abelian differential of the second kind with a double pole at the infinity w = oo. (So the 
above superpotential (4.38) has the sense of the Bloch dispersion law, i.e. the dependence 
of the energy w on the quasimomentum p.) The Hamiltonians of the KdV hierarchy can be 
obtained as coefficients of expansion of dp near the infinity. To obtain a complete family 
of conservation laws of the averaged hierarchy (3.1) one needs to extend the family of the 
KdV integrals by adding nonlocal functionals of u of the form 

i w'^dp, i w''-^dp, k= 1,2,.... (4.59) 

J tti Jbi 

As in (4.17) one can deform the above Frobenius structure on the moduli space M = 
Mg.no,-,nm by changing the differential dp, 

dp ^ dp = CA(f>t^ (4.60) 

for arbitrary constant coefficients. (The deformed Frobenius structure genericaly is well- 
defined only on a subset of M.) Particularly, if dp is a differential of the third kind on C 
then the "dimension" d of this model always equals 1. The corresponding hierarchy (3.1) 
is obtained by averaging a Toda-type system. 
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Here I consider the simplest example of such a deformation. Let us consider the 
3-dimensional family M of elliptic curves 



— 4{w — c)^ — g2{w — c) — — 4{w — c — ei){w — c — e2){w — c — 63) (4.61) 

with ordered roots ei, 62, 63. It is convenient to use the Weierstrass uniformization of 
(4.63) 

w = p{z) + c (4.62a) 

y = p'iz) (4.626) 

(I will use the standard notations [35] of the theory of elliptic functions). Let us use the 
holomorphic differential 

dp = (4.63) 

to construct a Frobenius structure on M (here p{u>) — ei). The corresponding Landau - 
Ginsburg superpotential is the Weierstrass function (4.62a) where one should substitute 
z = ajp/iri. The flat coordinates t^, t^, for the superpotential read 

f = -c+- (4.64a) 

UJ 

t^ = -l/uj (4.646) 

= 2'KiT where r = a;'/^, (4.64c) 

piuj') = rj — — p{z)dz. The charges of this manifold are qo — 0, qi — ^, q2 = d = 1. 

Remark. The above models with m = 0, g > can be obtained [39] in a semiclas- 
sical description of correlators of multimatrix models (at the tree-level approximation for 
small couplings they correspond to various self-similar solutions of the hierarchy (3.1)) as 
functions of the couplings after passing through a point of gradient catastrophe. The idea 
of such a description is originated in the theory of a dispersive analogue of shock waves 
[32]; see also [18]. 

More general algebraic-geometrical examples of solutions of WDVV were constructed 
in [44]. In these examples M is a moduli space of Riemann surfaces of genus g with a 
marked normalized Abelian differential of the second kind dw with poles at marked points 
and with fixed 6-periods 

^ Bi, i = l,...,g. 



hi 

For -Bj = one obtains the above Frobenius structures on Mg-no,...,nm- Unfortunately, for 
B ^0 the Frobenius structures of [44] does not admit a conformal invariance. 
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